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Abstr act 

An electromagnetic pulse (EMP) model Is developed which allows a 
quick computation of the time development of the electric fields qene- 
rated by a high altitude nuclear burst. The model Is based on the 
Karzas-latter high frequency approximation for high altitude EMP, which 
describes fields generated by Compton electrons interacting with the 
earth's magnetic field. The proper choice of a gamma time output func* 
tion, which can be integrated in closed form, and a small angle approxi- 
mation, made in the expressions for the Compton currents and air conduc- 
tivity, eliminate the time consuming numerical Integrations usually 
necessary in EMP models to compute the Compton currents and air con- 
ductivity. This results In a considerable savings In computation time. 
The mode' Is presented in a manner which is simple to use but still 
allows the vari ’■'on of the major theoretical parameters In the prob- 
lem. 

A. simplified model of electron collision frequency as a function 
of electric field strength is given which enables the model to predict 
accurate results for nuclear weapon gamma yields up to at least 100 Kt. 
The results predicted by this EMP model compare to within 5.5* with 
results from the Air Force Weapons Laboratory CHEMP computer code. 

The computation time using the presented model on a CDC &600 compu- 
ter Is typically 5 sec or less for a 5 shake computation period in 
steps of .1 shake. 

The mode! presented should be useful for both classroom instruc- 
tion and nuclear vulnerability/survivability studies and analysis prob- 
lems. 
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A CALCULATIOMAL MODEL FOR HIGH ALTITUDE EHP 

I . Introduction 

In recent years there has been a growing concern about the nuclear 
weapon effect known as the electromaonetlc pulse or the EMP. As the 
theory of EMP developed It was realized that there are actually several 
types of EMP which are characterized by the mechanisms which produce 
them. Types of EMP Include surface EMP, system generated EMP (SGEMP) , 
and high altitude EMP. A comprehensive discussion of the va 'ous types 
of EHP Is presented by Kinsley (Ref I). This report will address only 
one of these, that of high altitude EHP. 

The EMP Is basically a long range nuclear weapons effect wherein 
a high Intensity electromagnetic field Is radiated over a wide frequency 
band. The duration of this pulse Is typically on the order of shakes 
(1 shake » I0~ 8 sec). Specifically, high altitude EMP Is that produced 
by a nuclear detonation above 20 Km. 

The high altitude EMP problem is a subject of great Interest to 
the USAF due to its long range nature, and much effort goes Into pre- 
dicting the time development of the fields. These high Intensity EMP 
fields present a potential threat to weapons systems which rely on 
electrical or electronic components, such as communications, electronic 
counter measures (ECM), navigation, guidance, reconnaissance, and many 
others. The threat levels predicted by theoretical models have a di- 
rect Impact on the design of future Air Force weapons systems. These 
threat levels also help to predict the survivability of present weapons 
systems. 
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Th«> EMP generation models currently In use by the Air Force Wea- 
pons Laboratory (AFWL), are highly sophisticated and use large amounts 
of computer time. Other models available are very simple and give es- 
sentially order of magnitude results. The mode' presented in this re- 
port Is intended to provide an alternative betwetn these extremes, 
giving results close to those of the more advanced AFWL models but with 
very short expenditures of computer time. One important capability of 
this model is its ability to account for a preionization level due to a 
precursor burst. 

This model should be useful for quickly obtaining a meaningful 
estimate of the EMP environment for parameter variation or for use in a 
classroom situation to give a feel for the calculations involved In pre- 
dicting the generated EMP fields and showing the effects of a given 
parameter. 

The theory behind the high altitude EMP model presented is based 
on that of Karzas and Latter (Ref 2). Basically the Karzas-latter 
theory states that prompt gamma rays from the weapon produce Compton 
electrons within a specified region of the atmosphere known as the ab- 
sorption region. These Compton electrons are turned by the magnetic 
field of the earth to produce the radiated electromagnetic fields. 

There are several assumptions which are basic to the Karzas-latter 
model. These assumptions are: the earth's magnetic field is assumed 
to be uniform, the magnetic field lines have no curvature in the gamma 
ray absorption layer, the earth's surface Is assumed to be flat, the 
gamma rays are monoenerglc, each gamma Interacts with the atmosphere to 
produce only one forward directed Compton electron, the Compton elec- 
trons produced are also considered to be monoenergic and have a constant 
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velocity throughout their lifetime, and only the high frequency portion 
of the pulse is considered. 

One additional approximation Is made In this study. This is a 
small angle approximation in the trigonometric expressions for the Comp- 
ton currents and air conductivity.. This approximation, with the judi- 
cious choice of a function to represent the time dependence of the wea- 
pon yield, leads to closed form expressions for the Compton currents 
and air conductivity. This saves considerable computation time and 
greatly simplifies the calculations. For this approximation to be 
valid the cyclotron frequency u of the Compton electrons In the qeonag- 
netic field must be small enough so that the small angle approximations 
for sin ot and cos wt, where t Is in seconds, are valid. 

In all cases where quantities vary with the atmospheric density, 
an exponential atmosphere model is used. 

Factors not considered In the model are recombination of ions, 
avalanching or cascading of secondary electrons. X-ray effects, and 
self-consistent electromagnetic fields. 
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1 1 . Theory 


Overview 

Since the solution to the EMP problem is actually the solution of 
a classical electromagnetic theory problem, the derivation of the model 
equations reduces to putting Maxwell's equations into a convenient form. 
This is essentially accomplished by expressing Maxwell's equations in 
spherical coordinates and transforming to a retarded time frame. One 
must also develop expressions for the currents and conductivities of 
the system in the absorption region. The general equations describing 
the high altitude model of Karzas and Latter have been derived in great 
detail by Chapman (Ref 3). Only the major points of the derivation will 
be given here. The system origin is assumed to be at the burst point 
with detonation at time t = 0. This geometry is illustrated in Fig. I. 

The key points to be remembered in this model are: 

1. Each gamma ray gives rise to one downward traveling 
Compton electron. 

2. The electrons are turned by the earth's magnetic field 
giving rise to a centrifugal acceleration. 

3. The relativistic electrons radiate energy in their 
forward direction. 

A. The gamma rays and the EMP radiation travel at the 
same speed. This leads to constructive interference 
of the radiation from each of the electrons. 

P article Densities 

The gamma rays from >.ne nuclear weapon travel in a straight line 
to a point where they produce Compton electrons. At any given point r 


A 
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t’ic number of gammns which interact to produce Compton electrons Is 


expf-/,4C] 
v o * (r > 

£ " *orr 2 X (r) 


whore X(r) is the mean free path of gamma rays to produce Compton elec- 
trons, Y is the gamma yield of the weapon ir. electron volts (eV), and E 
is the mean gamma energy in eV, 

Equation (l) may also be called the radi3l distribution function, 

Y 

or an attenuation function for interacting gamma rays. The g- term is 
the total number of gamma rays available from the weapon. The 4«r 2 
term accounts for the divergence of the gamma rays as the radius r Is 
Increased while the remaining terms account for the reduction In gammas 
due to th-.ir absorption in the atmosphere, based on the mean free path. 

It is assumed that the gamma mean free path varies as the exponen- 
tial atmosphere. This gives the functional relationship between \ and 


Mr) * X g exp [(HOB - r cos A),-„ 


X q « gamma mean free path at standard pressure 
HOB « height of burst in Km above the earth's surface 
r ■ radial distance from the burst point to the point of Interest 
A « angle between the position vector r and the vertical 
S » atmospheric scale height 

With this assumption, Eq (1) can be Integrated and becomes 


, , Y 1 f S 

9 ( r ) - r T — TTTT exp i * , 7 

c liirr 2 X(r) ** Aq cos A 


r HOBn r f r cos A*\ , 1 1 

exp —J |_ exp — $ — ) - 1 J J 
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Mow if f(t) is the time distribution function of the weapon yield, 
the rate of Compton electrons, n c , produced at a given point r and time 


t is given by 


r-<« '0-0 


Each Compton electron produces through inelastic scattering events 
several secondary electrons which form the basis for the conductivity 
of the atmos- '•ere. As in the Kar 2 as-Latter approach, each Compton elec- 
tron is assumed to have a constant speed, V D , throughout the ranoe, R, 
of the electron which is a function of altitude. This allows the life- 
time to be expressed as R/V 0 . If each Compton electron produces secon- 
dary electrons at a constant rate, the rate of secondary electron, n s , 


production is 


dn s E c /33 «V 


where E c is the energy of the Compton electron and 33 eV is the average 
ionization energy per air molecule (Ref 2). 

Considering the differential current produced by the Compton elec- 
trons, it can be shown (Ref 3) that the Compton current and the number 
of Compton electrons are given by 


J c (t) = - e g(r) J ° V(t') f (x - t' + * — - *) dt' 

O C y 


n c (t) = g (r) / ° f (r - t' + dx' 


djfaifiia8aafewagfe.«.«8 
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i' “ the time since the creation of the Compton electron 
X(t") ■ the radial distance the Compton electron has traveled 
e « the magnitude of the electron charge 
The quantity t is generally known as retarded time. 

It then follows from Eqs (A) and (7) that the number of secondary 
electrons is 


r r° f G' - f ^ 


X(t")> . ..1 . . 

J dt ' I dt 


where q is E c /33 eV. 


Currents and Conductivity 

In the Karzas-Latter theory, the speed of the Co-ipton electrons Is 
considered to be a constant, however, there is an acceleration due to 
the geomagnetic field. The general equation of notion for an electron 
in this case is 


L. n yV = - e (E + V X B) - ra v c V 
dt 


"» = the electron rest mass 

V = the electron velocity 
E * the electric field 

B * the nagr.etlc field 
v c * the electron collision frequency 

Y “ (1 * (V 0 /c) 2 )*^ 
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If the relativistic motion of the electron is considered, only the 
V X B terra is iranortont and in sohcric.il coordinates, the expressions 
for the velocity components become (Ref 3) 


V r = V 0 (sin 2 6 cos ot + cos 2 0) 


V, -■ V„ (cos 9 sin 9 cos ot - sin 9 cos 9) (11) 

O O 


Vg = V 0 (sin 9 sin ot) 

where o i" the cyclotron frequency for an electron and is given by 


with 0 o the magnitude of the geomagnetic field. 
From Eq (10), X(t') is found to be 


X(t') = V 0 f sin 2 9 SI '"j U - T - + t' cos 2 bA 


The Compton currents may nov; be written as 


J c (t) - - cg(r) V 0 / ° (f(T)(cos 2 9 + sin 2 9 cos ot)l dt' (15) 


R/V. 

J C (t) * - eg (r) V / [f (T) sin 6 cos 8 (cos ut' - 0] dt' 


j c (t) » - eg(r) V / ° tf(T) sin 9 sin ot'] dt' 


9 
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i 

i 


where ft = — . 
c 

In a like manner £q (20) becomes 

T' - t' - (1 - 6) t" (77) 


An expression for the conductivity may also be found by using the 

equation of motion for the secondary electrons. These electrons are in 

the thermal regions with energies ranging from about 10-15 eV to the 

amoienc energy. It should be remembered for later use that the ambient 

energy of the secondary electrons is dependent on the electric field 

present. For consideration here, it is assumed that y \ I and also that 

the change of velocity with time is small compared to the other terms 

in Ea (9) so that nay be neglected. Also, with low velocities, the 

dt 

V X I term is small compared to the remaining terms and may also be 
neglected. Then the velocity of the secondary electrons is 

V „ - S_ E (28) 

mv c 

Using Eq (28), the current due to the secondary electrons is 

J s (x) - - e V n s (T> = n s (t) (29) 

Comparing Eq (29) to J s » oE, an expression for the conductivity 
is 

0 ( t ).sL„ s ( t ) ( 30 ) 

•' m 'e 

Equations (19), (23). (2*t), (25), (26). (27), and (30) provide the 
desired expressions for the Compton currents and the conductivity. 
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Fie ld Equot Ions 

Maxwell's equations is rational i zed MKS units are 


•* •+ aB 

V X E « - 22 . 

at 


•* + | ££ 
7 * B ■ V + ?fat 


» •* Qv 

7 .£ - — 

£ 0 


where q v is the total charge density and J is the total current density. 
In addition to these equations the continuity of charge requires that 


Combining these equations to separate E and E and transforming 

them into spherical coordinates and into the retarded tine frame (Ref 
-♦ •* 

3) the relations for E and S become 


1 ♦ ♦ 1 -*■ 

7 2 E + u r - — V.J + — Vq„ 

r CCg c Q 

3 7 2 1 3 , *, - . . n 

+ 57 .77 57 (rE) + u o (J * u r J r> j 


W-Po’*J*57 [ 77 57 (rB) 


+ 7- <V« - u 0 V J ’ 0 


In the Karzas-Latter model, only the time derivative portion of 
Eqs ( 36 ) and (37) are kept since the current variation with distance 
is slow compared to the variation in tine for the high frequency com- 
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ponents. Also the fields and currents vary r op idly in time. This 
approximation is valid for about 100 shakes. Th s sane hioh f reoucncy 
approx inat ion is used here. In addition, the radial component of the 
field is dropped since it is weak compared to the transverse components 
and contributes only a very low frequency sionol (Ref 2). The equa- 
tions for the transverse comnonents are 


It T f fIf (rE e.o ) + u o J «.f] 

(38) 


(39) 


(90) 

The currents in Eqs (33), (30), and (^0) are total currents. The 

total currents are given by 

J e,e “ J 0,2 + o(t) E 8,0 

(91) 

Substitution of Eq (91) into Eqs (33), (39), and (90) and 

tion over time gives 

intepra- 


I 7 57 (rEj + u 5 Jg + y 0 c(t) Eg = 0 

Iff? (r V + B O J 0 + **0 3(t > E fJ - 0 


I?! E s>'° 


Iff? (p ty + - J 9 + ~ o(t) E 6 ' 0 


m 

m 

m 


('.5) 
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Equations (f»?) and (^3) are In a form which can be solved. The 
ternt needed for solution of these equations arc the air conductivity 
and the transverse components of the Compton currents. The air con- 
ductivity nay be found by using Eqs (19) and (30)* The transverse com- 
oonents of the Conpton currents may be found by using Eqs (2*») and (25). 
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ill. MoJcl Cons Herat ions 


While the basic relationship expressions, have been developed be- 
tween E, J, and o, sene additional considerations are necessary in 
develooino the presented model. These include such items as the form 
of the function f(t), the average energies of the qamna rays and Comp- 
ton electrons, the mean free path of the na.mma rays and the range of 
the Compton electrons, and the collision frequency of the secondary 
electrons. The validity of the small angle approximation must also be 
examined. Each of these will be discussed separately in this section. 

When all of the items discussed in this section a-? combined with 
the equations given in the previous section, one obtains a straight- 
forward ard simple calculat ional model for predicting EMP environments 
from high altitude nuclear bursts. 



Time Dependence of Gamma Rays 

Many functions have been proposed to describe the time dependent 
output of a nuclear weapon, f ( t ) • These pulse shapes are often taken 
to be of the form (Ref *») 


f(t) = exp (- ?- - St) U(t) 


where U(t) is a unit step function havina a value zero for negative 
arguments and unity for ail others, and a and 6 are constants describ- 
ing the rise and fall of the pulse. Another form often used (Ref 5) I s 


(q + S) exp n (t - t p ) 

6 + a exp [ (a + S) (t - t 0 )] 
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where t 0 is the position of the peak value of the function. One fea- 
ture of these functions is that they cannot be intenrated : n closed 
form, and thus EMP codes which use the functions involve ime consuminq 
numerical integration schemes for calculating the Compt currents and 
the atmospheric conductivity. 

A basic idea in this investigation was to select functional forms 
of f(t) that could be integrated directly, thus yielding substantial 
savings in computation time. There are two constraints that the func- 
tional form of f(t) should follow. The first constraint is that the 
function be integrable in closed form. A second constraint is that the 
function should match reasonably well those currently in use for EMP 
calculations. 

Since each of the functions given in Eos {k&) and (A7) has an 
exponential rise and decay, prime consideration was given to functions 
having these same characteristics. With this in mind, the field was 
narrowed down to the following normalized functions: 

f(t) » lexp (-<*0 - exp (-St)] U(t) (A8) 

& - o 

... l [exp (St) - 1] U(t), t <. t 

f(t)= N 

[exp <Bt 0 ) - I] [exp - a(t - t„)] U(t - t 0 ) (1>9) 

where H is a normalization constant. 

It should be apparent at this point that the first constraint 
given could not be met without the simplification made to Eqs (16), 
(17). 08), and (20) by the small angle approximations In Eqs (21) and 
( 22 ). 
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After fui ther study, it was determined that Cq (43) would fit the 
functions qiven in Eqs (46) and (47) closely except for very narrow 

pulses. Figure 2 shews the shape of a typical long pulse. It is a plot 

of Cq (4£) with a = I0 7 and 6 <= 3.7 (10) 8 . The X narks indicate the 
values computed from Eq (46) with a = I0" 8 and 6 » 10 7 . Figure 3 shows 

a plot of a typical narrow pulse. In this case, Eq (A3) with a = 2 

(10) 8 and S = 3 (10) 8 , and Eq (46) with a « 2 (10)* 8 and 8 ■ 2 (I0) 8 

viere used. The values fron Eq (46) are again annotated by X marks. The 

curve computed from Eq (48) has been shifted in time away fron the ori- 
gin in order for the peaks of the two functions to match. This was done 

for an easier comparison. The result of this manipulation is to offset 

the times in the final calculations by the amount necessary to shift 
Eq (48). In both Figs. 2 and 3, the functions were normalised to a 
peak value of one f or comparison purposes. 

Equivalent comparisons using Eq (49) in place of Eq (48) indicated 
that Eq (43) would give a better overall versatility. Thus Eq (48) was 
chosen as the pulse shape fur tne sample calculations of this E.’IP model. 


Energies 

Two energy values are needed in the model, that of the gamma rays 
and that of the Compton electrons. Karzas and Latter (kef 2) suggest 
a I t'.cV gamma and a ! MeV Compton recoil energy. This implies no 
energy loss in the Compton creation process. Chapman (kef 3) uses 1.5 
MeV gamma rays and I.7.S HeV Compton electrons, which are the most ener- 
getic available from 1.5 HeV gamma rays. The values used in the AFWL 
models vary somewhat (Ref 4), but the most common value appears to be 
1.5 HeV gamma rays and .75 HeV Compton electrons. These Comoton elec- 
tron energies are actually the kinetic energy transferred to the elec- 
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tron from the gamma ray and assume the Compton electrons are only scat- 
tered forward. The encroy of .75 MeV from a 1.5 HeV qanma ray corre- 
sponds to the averanc encroy transfer in a Compton creation process 
(Ref 6). This Compton electron enerqy is qiven by 

T = (1 - a s /o c ) hv (50) 

where 

T = average recoil energy 

o s «■ average Compton scatterinq (non-absorption) cross section 
o c = average Compton collision (total interaction) cross section 
hv = energy of the gamma ray 

The values of a s and a c are from the Klein-liishina formulas (Ref 6). 

With these considerations, 1.5 HeV gamma rays are assumed while 
the Compton energy is left as a variable. 

A question directly related to the energies is that of gamma mean 
free path and Compton range. The mean free path is defined as the in- 
verse of the total cross section. For Compton creation the total macro- 
scopic cross section is the number of electrons times the absorption 
microscopic cross section. The absorption cross section is the differ- 
ence between the collision (total interaction) and scatterinq cross 
sections. Using this definition, the mean free path varies little for 
gamma energies between 1 and 1.5 HeV (Ref 6) so that the value sugges- 
ted by Karzas and Latter of 300 n at standard pressure Is a realistic 
value and is used in the model. 

Katz and Penfold (Ref 7) give relatlonshios between electron 
energy and range in any material. For electrons with energies less 
than 2.5 HeV the relationship is 
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where R is the ranqe in cm, E 0 is the energy in MeV of the electron, 
and 0 is the density in mg/cm 3 of the material of interest. For a .75 
MeV electron and an air density of 1.293 mq/cm 3 at standard pressure, 
the range is 2.23 n. For a 1 MeV electron the range is 3-18 m, which 
shows an increase of almost one meter. Since the Compton recoil energy 
is taken as a variable and the range variation is considerable, Eq (51) 
with an air density of 1.293 mg/cm 3 is used in the model to compute the 
range of the Compton electrons at standard pressure. 

Another factor which is a function of Compton electron energy is 
the velocity of the Compton electrons. The. energy used for the Compton 
electron is considered to be the kinetic energy of the electron. The 
total energy is the sum of the kinetic and rest energies. The relati- 
vistic relation here is 


1 


E + ro 0 » m 0 y (52) 

where E is the kinetic energy in eV and m 0 is the rest energy in eV. 
Equation (52) may be solved for the Compton electron velocity to get 


Vo 


c 



m o 

E + m D 



(53) 


Equation (53) is also used in the model. 

The energy from the prompt garma rays is assumed to be deposited 
in an absorption region from 20 Km to 50 Km above the earth's surface. 
This is based on calculations by Latter and lelevier (Ref 8). 
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Compto n El ect ron Lifet ime 

The lifetime of a Compton electron, as previously defined, is 
given by the quotient of the electron ranoe P. and its constant speed 
V 0 . The ranoe at any given altitude is found by exoonent ia I ly scaling 
the standard pressure value found from Eq (51) and the sea level air 
density. Since the high frequency approximation is only valid for 
about 100 shakes, the lifetime should not be allowed to exceed this 
length cf time. 

Small Ann i e Approximation 

One inportant factor for this model is the validity of the small 
angle approximations leading to Eqs (2f» ) , (25). (26), and (27). This 
can be somewhat arbitrary, depending on the amount of error one Is 
willing to accept. For the use of this analysis, a 1% error will be 
considered acceptable. With this limit the maximum value for the solu- 

i 1 

tion of sin A ■ A is A = .36. i i 

The maximum geomagnetic field in the qar.ma ray absorption region " 

is usually taken to be a .6 gauss (Ref 4). A minimum kinetic energy ; r 

J ; 

of .5 MeV will be assumed for the Compton electrons. The time which . . 

satisfies : 's 

or = .yh (54) 

is the maximum time duration over which the small angle approximation 

is val id. ; 

Equation (13) gives the expression for j. This will be a ma :imum 
when B 0 is a maximum and y is a minimum. For y to be a minimum, the 
particie speed must be a minimum. The speed of a .5 MeV electron is j 

I v 

.863c. With these factors u is 5.33 (I0) 6 . Solving for t in Eq (5b), | : 
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the result is t ° 6.*i shakes. This is however a very worst case. In 
(rest situations the electron energy will be considerably greater than 
•5 MeV. For on electron of .75 McV with a speed of .9l**2c, the time 
limit is t s 7.9 shakes. If the geomagnetic field Is .3 gauss, which 
is often the case, the time limits above will be doubled. 

The minimum length of the confidence interval should be C.k 
shakes. As the electron energy increases, the length of the confidence 
interval will also increase. 


Electron Collision Frequenc y 

In previous EH° models two general approaches have been taken to 
find a value for the air conductivity. The first uses Eq (30) and 
assumes a constant collision frequency scaled to an exponential atmo- 
sphere. This approach tends to limit the usefulness to a rather nar- 
row range of gamma yields if accurate answers are desired. This basic 
approach is the one suggested by Karzas and Latter. 

The second approach is the use of the electron transport equations 
(Ref 9) to solve for the secondary electron velocity. This velocity is 
then used to compute the secondary currents directly. These equations 
however, require an iterative solution which is very time consuming. 

Since the model developed here is desired to be as simple as 
possible to use and understand, but still give good results over as 
wide a range of gamma yields as possible, an intermed iate approach Is 
desirable. 

When a secondary electron is created, it has an energy somewhat 
higher than Its ambient energy. The secondary electrons transfer this 
additional energy to air molecules through both elastic and inelastic 
collisions. The time required for this process to take place Is known 
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os tic thermal i zat ion time. Eaun (Ref 10) has mothered data to show 
that ti>e thcrr.al i rat ion time for low energy (otout 10 cV or less) elec- 
trons at sea level is approximately a nanosecond. Usino an exoonential 
atnosphere model, this thermal izat ion time corrcsnonds to approximately 
2 shakes at 20 Km and to over a microsecond at 50 Km. The thermal iza- 
tion time in the 20 Km case is on the same order as the time scale of 
interest (a few shakes) for the high altitude EM D problem. The scale 
height used for the exponential atmosphere is taken as 7 <m. This value 
is based on a curve fit of data from the U.S. Standard Atrosphere of 
1962 (Re r II). This scale height of 7 Km is used throuohout the model. 

Baum also shows the effect of an applied electric field on the 
ambient electron energy. The secondary electron ambient energy in- 
creases with increasing field strengths above I0 1 * vo'ts/r.. For elec- 
tric field strengths below I0 h volts/m there is little or no effect. 

If the thermal ization process is assumed to be a linear function 
of tine and the lower energy linit of the secondary electrons is as- 
sumed to be a linear function of the electric field strength, these 
effects can be easily incorporated into the E U .P rodel beinn developed. 

It is however, necessary to determine the coefficients of the assumed 
functions. 

The tine dependent function is based on the thermal izat ion tine 
at 20 Kn. The coefficients for the electric field dependent function 
are based on an empirical fit to the peak electric field values from 
the AFWL CHEHP model (Ref E). In order to get an acccntable fit, It 
was found that the field dependent function had to be split into two 
linear functions. Since the electron energy and collision frequency 
are directly related (Ref 9), the linear relations described above may 
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be expressed directly as electron collision frequency. The functions 
described arc 


-2.5 do) 20 T + k.5 (I0) 12 


(55a) 

0.k3 (10) s E + 1.6 (I0) 12 

E < s do) 1 * 

(55b) 

6 (10) 7 E + 0.6 (I0) 12 

E > 5 (10)" 

(55c) 


where t is retarded tine in seconds and E is the electric field 
strength in volts per neter. 

The collision freouency at sea level is taken as the maximum of 
Eos (55) and 2.S ( 1 0) 1 2 . In no case is the collision frequency allowed 

to exceed k.k {1 0) 1 2 . The collision freouency is then scaled using an 
exponential atmosohere. The lower limit value of 2.3 (10) 12 was deter- 
mined empirically while the upocr limit value of k.k (lo) 12 results 
from the collision frequency being essentially independent of the elec- 
tric field for extremely high fields (Ref 2). 

P reion i zat ion 

The problem of preionization by a precursor burst may be handled 
quite easily by the model developed here. This is because of the 
choice of the oamna yield function of Eq (kC) and the small angle ap- 
proximations of Eqs (21) and (22), which allows Eq (I?) to be Inteqra- 
ted in closed form to get the nunber of secondary electrons as a func- 
tion of retarded time for any altitude. The total air conductivity 
for a burst with a precursor burst is comouted from the total number 
of secondaries from both bursts at the time of interest. The colli- 
sion frequency should be based on the time since the precursor burst 
and the electric field strength due to the main burst at the time of 


a 


i 



interest. 
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it is In the calculation of the number of secondary electrons from 

a precursor burst where the tine confidence limit of the snail angle 

approximations is greatly exceeded. The term of interest here is 

LI — , — in Eq (20). Since this term and its small amle approximation 

u 

t" both have maximum on the order of I0" 6 , and recombination Is a 
noticeable factor, on preionization time scales, any error introduced 
by the small angle approximation should be tolerable. 


Calculation Method 


The time development of the EKP signal may be found by the numeri- 
cal integration of Eqs (A2) and (A3) over the absorption region along 
the line of sight between the burst point and the target for retarded 
times from 0 to i. The distance from the burst point to the top of 
the absorption region along the line of sight is called RMIti while the 
corresponding distance to the bottom of the absorption region is called 
RMAX. The distance to the target from the burst point is r. It should 
be noted that above and below the absorption region there are no cur- 
rents or air conductivity so that the only contribution to the electric 
field is over the absorption region. The integration limits are then 
RM1N and RMAX as shown in Fig. I on page 5. RMAX will equal r If the 
target is in the absorption region. 

Because a numerical integration Is necessary, this model is best 
run on a computer. The results given in this report were calculated 
through a computer code usino the concepts and methods detailed here 
applied to the theory given previously. The computer code calculation 
will be referred to as HAEMP in the next section. 

For points below the absorption region, where both the Compton 


currents and air conductivity are zero, Eqs (A2) and (A3) have soiu- 
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hold, the first term in Eqs (42) and (43) may be neglected and they 


Eg » - Jg/o(r) (61 

Eg»-J§/o(t) (62 

the values found from Eqs (6l) and (62) may then be used with Eq (53) 
to find the electric field at the target. 
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hold, the first term in Eqs (42) and (43) may be neglected and they 


Eg » - Jg/o(r) (61 

Eg»-J§/o(t) (62 

the values found from Eqs (6l) and (62) may then be used with Eq (53) 
to find the electric field at the target. 
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IV, P esul ts 

To test the validity of the liAEMP model, the ccmnuted values of 
this model were compared to available equivalent model values from the 
AFWL code for EMP calculations t.nowr as CIIEttP (P.ef A). 

The basic set of conditions used for these calculations was: 
target location = nround zero 

height of burst = 100 Km 

geomagnetic field = .3 nauss or 3 ( « 0) ~ s wb/m 2 

inclination angle = 0.0 degree 

Comptcn electron recoil energy = .75 HeV 

The gamma yield was va : ed and the number of steps taken for the 
numerical integration was varied according to the gamma yield, with 
more steps taken for the higher yields. Other parameters varied for 
examination of ocak field values were burst height, geomagnetic field, 
and pulse shape. A preioni zat ion level was also considered. 

The available data from the CHEMP(ll) code, which is CHEMP run with 
non-self-consistent calculations, was computed using a pulse of the 
form of Eg (46) with a » 10* 8 and S = 1C 7 and the same geometry as 
given above. This pulse shaoe is almost identical to that of Eq (AS) 
with a 10 7 and 6 = 3.7 (10) 8 . See Fig. 2 on page !£. Using this 
pulse shape, a range of gamma yields from .01 Kt to 100 Kt was used to 
calculate the EMP field values. The peak field values are plotted in 
Fiq. 4. The values taken from the CHEMP(N) code are annotated by X 
marks. These peak field values shew a maximum difference of 5.5% 
around the I Kt case, and a difference of less than 2% for all other 
known cases. 
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The complete pulse out to i < ihak.es for the .29 Kt case is shown In 
rjq. 5. '■.gain, values from the '.HW ('0 code are indicated by X marks. 
There is excellent agreement for the first two shakes after wh'ch time 
the IIAEKP values fall off more rapidly. This more rapid fall is most 
likely caused by two factors. The first is not considering recombina- 
tion of the secondary electrons, which would tend to reduce the con- 
ductivity and allow the fields to faP off more slowly. The second is 
the slew breakdown of the small angle approximation in comnuting the 
Compton currents. 

The height of burst was varied to observe the effect on the peak 
fields for the various yields. This data is plotted in fig. 6 for 
altitudes of 60 Krt, 100 Km, 200 Km, and 300 Km. The neneral trend of 
Fig. 6 indicates that the field strength increases as the height of 
burst increases only for the larger yields, for the lower yields, the 
field strength decreases as the burst height increases. Examination 
of the hypothetical curves in fig. 6 implies that the maximjm EMP sig- 
nal generated by any given yield depends on the height of burst, for 
example, using the curves in Fig. 6, it could be concluded that the 
maximum EHP from a I Kt gamma yield would occur for a height of burs 
somewhere between 60 Km and 200 Kn and for a . 1 Kt oamma yield, the 
burst height would be less than 100 Km. 

A test was also run to determine the effect of a narrower gemma 
output pulse. A pulse used by the AFVfl in the CHEMP code for this 
ourpose is of the form of Eq (A6) with a - Z (10) 8 and 3 * 2(10) 8 . 

An essentially identical pulse can be achieved with Eq (W) using 
a * 2 (10) 9 and 3 = 3 (10) 8 . See Fig. 3 on page 19- A plot of peak 
fields vs yield Is shown in Fig. 7 for both the narrow pulse and the 
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The complete pulse out to i < ihak.es for the .29 Kt case is shown In 
rjq. 5. '■.gain, values from the '.HW ('0 code are indicated by X marks. 
There is excellent agreement for the first two shakes after wh'ch time 
the IIAEKP values fall off more rapidly. This more rapid fall is most 
likely caused by two factors. The first is not considering recombina- 
tion of the secondary electrons, which would tend to reduce the con- 
ductivity and allow the fields to faP off more slowly. The second is 
the slew breakdown of the small angle approximation in comnuting the 
Compton currents. 

The height of burst was varied to observe the effect on the peak 
fields for the various yields. This data is plotted in fig. 6 for 
altitudes of 60 Krt, 100 Km, 200 Km, and 300 Km. The neneral trend of 
Fig. 6 indicates that the field strength increases as the height of 
burst increases only for the larger yields, for the lower yields, the 
field strength decreases as the burst height increases. Examination 
of the hypothetical curves in fig. 6 implies that the maximjm EMP sig- 
nal generated by any given yield depends on the height of burst, for 
example, using the curves in Fig. 6, it could be concluded that the 
maximum EHP from a I Kt gamma yield would occur for a height of burs 
somewhere between 60 Km and 200 Kn and for a . 1 Kt oamma yield, the 
burst height would be less than 100 Km. 

A test was also run to determine the effect of a narrower gemma 
output pulse. A pulse used by the AFVfl in the CHEMP code for this 
ourpose is of the form of Eq (A6) with a - Z (10) 8 and 3 * 2(10) 8 . 

An essentially identical pulse can be achieved with Eq (W) using 
a * 2 (10) 9 and 3 = 3 (10) 8 . See Fig. 3 on page 19- A plot of peak 
fields vs yield Is shown in Fig. 7 for both the narrow pulse and the 
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Fig. 6. Peak Electric Field as a Function of Curst Height 
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original long pulse. This data points to n saturated peak field 
strcnath of A. 32 (10) 1 * volts/m for the 100 Yn burst case in a .3 oauss 
geomagnetic field. The narrow pulse curve shows that the peak fields 
ore higher for all yields. The data from Fin. 6 however shows that 
the saturated field value preoictcd in Fin. 7 does not hold for all 
burst heights. 

Another factor which should have a noticeable effect on the EMP 
generated is the presence of a precursor burst. For this test a .03 
Kt qamna yield burst was assumed to have bpen detonated 10 nsec prior 
to the main burst. The pulse shape used to describe the qamma output 
is the initial long pulse ana Is used for both bursts. The peak field 
results are plotted ir, Fie, G. The results of this test were au.te 
dramatic. This data shows that the pear, field of 1.6 (ID) 4 volts/n 
from the .03 Kt precursor burst is not equaled until the main burst has 
a yield of nearly 1 Kt. In general, the neak field values with this 
level of oreicni7ation are tho same as If the weapon yield was one 
tenth of its actual yield. As the yield of the main burst goes up, the 
field reduction becomes less. At a yield of .1 Kt the field is reduced 
by over a factor of A. At .25 Kt the factor drops to 3 and to 2 for 
the I Kt case. When the yield finally reaches 100 Kt, the reduction 
factor has dropped to 1.1. 

The results with preionization compare favorably with available 
results from the AFWL (Ref A). However, since a slightly different 
geometry and some additional source modeling were used, a direct com- 
parison cannot be made. 

A last major factor to consider for this model is the geomagnetic 
field strength. The peak field values as a function of yield are 
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plotted in Fig. 3 for the .3 gauss ana the .6 gauss ocenannctic fields. 
The larger geomagnetic field resulted in apnroxinately a tines 2 In- 
crease in the field strengths for all vields considered. 

One renaining result should be mentioned. This concerns the tine 
of the peak field value. Since for these calculations the fields were 
computed in steps of .1 shake, the peal field tine will he accurate 
only to this Unit. The peak for the .25 Kt case shown in Fig. 5 oc- 
curs at 1.7 shakes for the IIAEHP node! and at an estinated 1.3 shakes 
for the CHEMP(N) code. As the yield goes up, the peak field value 
occurs earlier. For yields above .25 Kt, the peak occurs between 1.1 
and 1.3 shakes for all cases examined. For all yields and altitudes 
considered, the peak always occurred within 8 shakes, with the majority 
falling between 1 and 2 shakes. 

Since the model presented is supposed to be a quick computational 
tool, some estimates on the calculation tines that can be cxcccted are 
appropriate. The tines given here are the average tines used in com- 
puting the results given above and ore based on a computer code run on 
a CDC 6600 computer. The average computation tine was 1.5 sec for a 
5 shake calculation of the electric field development computed every 
.1 shake using 50 steps in the numerical integration. Increasing 
either the parameter of 50 steps or 5 shakes by any factor increased 
the computation tine by the sane factor. The presence of a precursor 
burst had no noticeable effect on computation tine. When the approxi- 
mations of Eqs (61) and (62) were valid, the tine dropped to less than 
.2 sec for a 5 shake calculation. 

An interesting observation related to the numerical integration is 
that as Eqs (61) and (62) become valid, the numerical integration could 
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V . Conc lusi o ns and RecovroncJo t i ons 

This report has presented a ealeu lot ional model f or the hioh al- 
titude EtfP. This model falls between the very simple and niqhly sophis- 
ticated models currently available. The fnodel is easy to understand 
but allows the variation of the major theoretical parameters. The model 
is also simple enough to use so that a preionization level car be in- 
corporated with a very minimum of effort. The results given by this 
model are in excellent agreement with the available data from the very 
sophisticated models usina the same geometry. The computation tine is 
on the order of a few seconds. 

Uses 

There are two primary uses seen for this model. The first Is for 
sensitivity determinations for a given parameter, such as precursor 
yield in preionization studies, height of burst for a fixed yield, or 
changes in the magnitude and inclination of the neonagnetic field. 
Calculations using the HAEMP model can help to evaluate the relative 
importance of the many variables encountered in EMP problems. It is 
useful in this role because it oives quick answers which arc reasonably 
accurate. 

The second use is as a classroom model. The basic physics of the 
high altitude problem is present in an easy to use form which demon- 
strates how the currents, conductivity, and electric field interact 
with each other. Because the calc si t ions reauired to get an answer are 
minimal and straightforward, the physical processes are not lost in the 
search for a result, as can happen all too easily with hiqhly sophis- 
ticated models. Since the results are in relatively good agreement with 
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the sophisticated models, the student can leave the classroom with a 


oood concept of the overall problem. 


Limitations 


The model does have its limitations, however. These are due to the 
various assumptions and approximations used to arrive at the final 
calculat ional model. These assumptions and approximations should be 
kept in mind when using the model. If this is done, the results from 
any problem can be reasonably interpreted. 

Two of the error sources are the small angle approximations and the 
assumption of a flat earth. For the small angle case the error Is 
nealigible at very early times but increases as the time increases. 

In the flat earth assumption, any error is relatively constant for a 
particular geometry, but changes when the geometry changes. This error 
increases as the ground inane of the target moves on the earth's sur- 
face away from the ground zero ooint. The other assumptions also add 
some error to the results, but if the effects of these assumptions are 
renembered, the KAEMP model can be helpful for the solutions of a wide 
range of EHP problems. 


P.ec omnenda 1 1 ons 

There are several areas of investigation which might lead to Improve- 
ments in the HAEMP model. One of the more important possible investi- 
gations is to determine if the electron transport equations used In the 
AFWL CilEHP model to describe the secondary electrons can be greatly 
simplified and used In the HAEIIP model. Another area is the possibility 
of allowing recombination and cascading of secondary electrons. The 
use of self-consistent electromagnetic fields Is also a possible source 


1(1 


of Improvement. This change however would have to he examined very 
closely to determine the effects it would have on the validity of the 
small anqle approximations, which ore vital to the IIAEMP model. An 
additional possibility is the consideration of small anole scattering 
of the Compton electrons. This could possibly be included by a random 
elimination of some Compton electrons for Compton current production 
while still allowing all Compton electrons to produce secondary elec- 
trons. 

Each possible change to the HAEHP model should be examined with 
consideration to the tradeoffs between more complete solutions and the 
continued ease cf use and understanding of the model. 
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HApl P ^ odp I Stun nary 


The user of the HA£MP mode! must first determine c he following 


parameters in MKS units: 


the gamma yield Y of the weapon in MeV. 


the height of burst (hOB) of the weapon. 


the magnitude 3 0 and angle of inclination of the 


geomagnetic field in the absorption region. 


the target location in spherical coordinates with the 


origin at the burst point and wi th trc polar axis 


parallel to the geomagnetic field lines. 


the kinetic energy E c of the Comotor, electrons In MeV. 


the parameters defining the gamma output pulse shape. 


the atmospheric scale height S to be used. 


From these parameters, the values for the annle A and the two dis- 


tances ®?1IM and RMAX used for integration over the absorption reoion 


can be found from geometry. The angle from the vertical to the line of 


sight from the burst point to the target is the angle A, as shown in 


-ig. 1 on page 5. The value of RHIM is the radial distance from the 


ourst point to the top of the absorption layer (50 kn above the earth’s 


surface) along the line of sight to the target. The value of RMAX is 


the corresponding distance to the bottom of the absorption layer (20 


Km above the earth’s surface). See Fig. 1 on pane 5. 


The electric field development from tines 0 to t in steps of At in 


the retarded tine frame is found by the numerical integration of 


1 7 If (rE e> + u o J 8 + "o o(t) e 8 " 0 
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"fir (rE 0> + "c \ + »c «< 1 > E„ • 0 (* 3) 

over r from RMIM to RMAX for each tire increment. The initial condi- 
tions on these equations are that E 9 a rd E ? are zero at t = 0. The 
number of steps taken over the inteeration interval will depend on the 
yield of the weaoon and the accuracy desired, ho less than 50 steps 
should be used or more than 500 steps necessary for gamma yields up 
to 100 Kt. 

The solutions to Eqs (k2) and (A3 ) will give the electric field 
components at r <* RMAX. The total field strenoth Is oiven by 

E = [(E 0 ) 2 + (Eg) 2 ]^ (63) 

and the field strength at the target is given by 

, RMAX E 

E (target) « (53) 

'target 

The pe.iK electric field value over the time interval 0 to r Is the 
maximum value found using Eq (53). 

In the special case where the conditions 

If •!? (rE eyi « “0 i J e,jii (59) 

cf if? (rE 8,e>l c< b 0 o(t) E 9>0 (60) 

hold, the electric field corcoonents are alven by 

E fl - - Jfl.Mt) (61) 


Eg ‘ • Jjj/d(*) 


dl 
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fhc values from tqs (61) and (6?) nay then be used with Eqs (63) 
ar.d (53) to find the electric field at the toroet. 

In order to integrate tqs (A?) ar.d (A3) or to compute tqs (61) 
and (62), the terms J§, J§, and o(t) are needed as functions of time i 
and distance r. These are given by 


j!i = eq(r) V c sin 0 cos 8 — j t' j f(T) dt' 


Jg = - og(r) V Q sin 8 a | t' f(T) dt' 
o 


o(t) = n^ n S< T > 


CXD * - T 

£ c lerr 2 X(r) f X 0 


, HOB ^ ^ r cos A v . 

exp - ' — ) |_exr, ^ — T~ j ' 1 i } 


f(T) = [exp <-aT) - exn (-bT)] U(T) 

b - a 

Q V_ t R/V 

n «(t) ■ — — g(r) j [/ f(T') dT"] dT- 


t = t - U-e) t' 


T' = t' - (1-S) T' 


withe - ar.d q = E c /33 eV 


The Compton electron cyclotron frequency u is given by 


03 ) 


0£P/PH/7 r ,-l3 


r 


F: 




! 

i 


i 



1 


L 


The Conoton electron velocity V 0 is oiven by 


V 


o 


c 


( — !V 

^E c + m 0 > J 


The Comotoo electron range at sea level is given by 


( 53 ) 


R = (612 £ c <’- 26 5 - ’W' E c ) )/D (5 , , 

where 9 = 1.293 mq/cm 3 . This range is then converted to meters and 
scaled to the appropriate value for the point r by multiplying the com- 
puted value by 

exp {(HOB - r cos A)/S] (66) 

Then the Compton lifetime P/V 0 is oiven by the quotient of the correc- 
ted P and V 0 , but is never allowed to exceed one usee. 

The value of A(r) is given by the product of (300 m) and Eq 

(63). 

The collision frequency v c at sea level is given by the maximum 
of 


v c = -2.5 do) 20 t + 6. 5 do) 12 


(55a) 

v c = 0.63 (i0) 8 E + 1.6 (10) 12 

E < 5 do) 1 - 

(55b) 

v c = 6 (10) 7 E + 0.8 (10) 12 

E > 5 (10)“ 

(55c) 


or 2.3 (10) 12 where t is the retarded time in seconds and E is the elec- 
tric field strength in volts per meter. The collision frequency is not 
allowed to be greater than 6.6 (1 0) 12 . The value of the collision fre- 
quency is recomputed after each tine step in the development of the 

<l3 




,J 


I 
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electric field. The sea level value is also scaled for altitude by 
dividin'} by Eq (of* ) . 

The intcqrals in Eqs (24), (25). and (19) can be evaluated in 
closed form. The steo function II (T) in Eq (43) results In two expres- 
sions for each of the inteorals. The first set of exnrcssions Is 
used when r — (? - 2) (R/V 0 ) <.0 and the second when x * (I - S) (R/V Q ) 

> 0 . 
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J c (t) => eo(r) V sin 6 cos 0 


(I - K ) 3 


{exp (-at) {exp (aL) [(aL) 2 - 2aL + 2) - 2) 


exp (-bt) (exp (bL) [ (bL) 2 - 2bL + 2] ~r) 

b* 


J§(t) = - cp(r) V 0 sin 6 w 


b - a (, . B )2 


{exp (-at) {exp (aL) [aL - l) + i] 


- exo (-bt) (exp (bL) f bL - 1) + 1] r) 


o(t) = -~g(r) - rr n,r s 


(|- [at - 1 + exp (-aL) - [exp (aL) - l][exp (-at) - exp (-aL)]] 


- bL - 1 + exp (-bL) - [exo (bL) - l][exp (-bT) 
b 


- exo 


where L * (1 - 3)(R/V 0 ). 
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